Estimate of the long-distance contribution 
through b ^ sip to the Bg 77 decay rate 



G. Hiller*, 

Deutsches Elektronen-Synchrotron DESY, Hamburg 



E. O. Iltan t 

Physics Department, Middle East Technical University 

Ankara, Turkey 



Abstract 

The Bg — > decay is modeled through the inclusive h ^ sip decay. Using the Vector 
Meson Dominance model, the amplitude for the chain process — > (pip — > (/)7 — > 77 is 
estimated and it is found to be at most 4% of the corresponding amplitude from the O7 
type LD contribution. The intermediate amplitude for the process Bg 07 is compared 
with the corresponding one obtained by a different approach based on the interaction of 
the virtual charm quark loop with soft gluons Q. Both amplitudes are found to agree 
within 10%. Further the influence on the branching ratio B{Bs — > 'yj)sD+LDo7 from 
inclusive b ^ sip is estimated as less than 1%. 
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1 Introduction 



Rare B-meson decays are one of the main research fields in recent particle physics. Since 
they occur only at loop level in common fiavour changing neutral current (FCNC)- forbidden 
models, the theoretical and experimental investigations provide precise tests of the Standard 
Model (SM) and possible new physics beyond. An interesting candidate for rare decays is 
Bg — s> 77, which has a rich final state since the two photon system can be in a CP-odd or a CP- 
even state. In the literature, 77 has been studied in the lowest order in the framework 
of the constituent quark model and found to consist in the effective Hamiltonian theory 
which is obtained by integrating out heavy particles (top quark and W boson in the SM) 
from the full theory, of contributions from O7 and four-quark operators. The calculation for 
Bs 77 decay with leading logarithmic QCD corrections including O7 type long distance (LD) 
effects and in an approach based on heavy quark effective theory to model the bound state was 
recently done in 0. 

In the present work we continue and estimate the additional LD effect due to the dominant 
four-quark operators Oi and O2 (see eq. (0)) through the B^ — > (pip — > ^7 — > 77 chain decay. 
We use at quark level b ^ sip followed by the b —>■ sj decay and we pass to the hadronic 
level using the transition form factor -Fi(O) from the amplitude A{Bs — 0) 0,10. For both 
the conversions ip ^ 'j and ^ 7 we impose the Vector Meson Dominance (VMD) model [^]. 
The conversion ip ^ 'y needs further manipulation because of the strong contribution of the 
longitudinal part of the tp meson. We extract the transverse part using the Golowich-Pakvasa 
procedure Q, f^. Further we calculate the Oi,2 type LD effect to the Bs 07 decay using 
the method given in , namely, by taking into account the virtual c-quark loop instead of the 
hadronization of the cc pair. This procedure was originally applied to estimate the LD effect 
in B ^ K*'y decay and uses operator product expansion and QCD sum rule techniques. We 
compare the amplitudes obtained by these two different methods and show that they are in 
good agreement within the errors of the calculation. Finally we estimate the influence of the 
Oi,2 type LD contribution on the branching ratio for Bg —>■ 77 decay. 
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The paper is organized as follows: In section 2, we calculate the Oi 2 type LD contribution 
to Bg — ^ 77 due to the chain process Bg — > (pip — 07 — 77 using the Gordon decomposition 
and further the VMD model. We compare our amplitude for the decay Bg (pip (p'j with 
the Oi,2 type one for Bg — ^ 07, calculated by virtual charm loop-soft gluon interaction [^]. 
Section 3 contains a discussion of the branching ratio B{Bs — > 1i)sd+ldoj+ldo2 "^^ere the 
subscripts SD, LDpr, LD02 denote QCD short distance (SD) contributions and LD ones due 
to Bs ^ (p'~i 77 and Bs — ^ (pip ^ 07 — *• 77 decays, respectively. 

2 The chain process Bg (pip ^ (/)7 ^ 77 

We first consider the additional contribution to 6 — 57 from b — > sipi 57, where ipi are all cc 
J = 1 bound states, see fig. |l|. The relevant part of the effective Hamiltonian describing this 
process is given as 

neff = A^v:MCiif^)OM + C2(/i)02(/i)), (i) 

with the dominant four-Fermi operators 

1 — 75 _ 1 — 75 

01 = Sa7M 2 ^'^ ^'^^^ 2 ' 

^ _ 1 - 75 - 1 - 75 , .^s 

02 = s^^^^c c^^^^b . (2) 

Here a, (3 are SU{?)) colour indices and V^*^ are the relevant elements of the quark mixing matrix. 
The initial values of the corresponding Wilson coefficients are Ci(mvy) = and C2{Tny/) = 1. To 
include leading logarithmic QCD corrections we evaluate Ci,2(/^) at the relevant scale, fi ^ rrib 
for i?-decays, and this takes into account short distance effects from single gluon exchange. 
The analytical expressions can be found in 

Further we have used the unitarity of the CKM matrix V*^Vcb = —V^*Vtb — V*gVub and have 
neglected the contribution due to an internal u-quark, since V*gVub « Vts^tb = ^t- 

Using factorization, we obtain the inclusive decay amplitude for the process 6 — >^ s?/' as 

^(6 ^ siPih, e"^)) = -iCU{ml)m^-sr{l - l,)hel . (3) 
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Figure 1: The diagram contributing to Bg (j)^ 07- 

Here 

C = -^XtaM (4) 

with, assuming naive factorization, 

«2(/.) = CM + ^, (5) 

where Nq = 3 in colour SU{3) and ki, are the momentum and the polarization vector of the 
ip, respectively. In eq. (|^) we used the matrix element 

< 0\c'y,c\ij{h, e'^) >= /v,K)m^eJ . (6) 

At this stage there is a critical remark about factorization in order, concerning the value of 
02 (/i) used. The decay under consideration is a class II decay following the classification of [§]. 
In general eq. is written as 

a," = + ^) [1 + eM] + CMe2ifi) , (7) 

where ei(/i) and e2(/i) parametrize the non-factorizable contributions to the hadronic matrix 
elements. a2^^ takes into account all contributions of the matrix elements in contrast to 0.2 (/u), 
which assumes naive factorization ei(/x) = e2(/i) = 0. Especially e2(/i), which is the colour 
octet piece, has sizable contributions to naive factorization in class II decays Furthermore, 
the additional problem is not to know the correct factorization scale. In order to include the 



non-factorizable corrections we use in the definition of C in eq. @ the effective coefficient , 
which is determined experimentally from the world average branching ratio of i? ^ K^*^ip as 

i 

a'J^ = 0.21 . (8) 

This choice restores the correct scale and is n independent. Writing the Ansatz 

4// = Ci(m,)+eC2(m,) , (9) 

it follows that C, ~ 0.41 with Ci{mh) = —0.25 and C2{mh) = 1.11 for the input values given in 
Table |1|. For comparison, naive factorization would give a2{mb) = 0.12. 



Parameter 


Value 


rUb 


4.8 (GeV) 


rric 


1.4 (GeV) 


"em 


129 


Xt 


0.04 


rrit 


175 (GeV) 


mw 


80.26 (GeV) 


mz 


91.19 (GeV) 


a(5) 


0.214 (GeV) 


as{mz) 


0.117 



Table 1: Values of the input parameters used in the numerical calculations. 

Now our aim is to replace the ip meson with the photon 7 and to construct a gauge invariant 
amplitude. This can be done by killing the longitudinal component of the meson so, that 
is changed to the polarization vector of the photon 7. For this we use the Golowich-Pakvasa 
procedure with the help of the Gordon identity, namely 7^7^ = g^a — 'i(^tj,a- We start with 
the vertex 57^(1—75)6 and using the equation of motion = rribb and momentum conservation 
p = p' + ki, we get 

S7^(l - 75)6 = — {s7^ /(I + 75)^ + s7m A(1 + 75)^} , (10) 
rUb 

where p,p' are the momenta of the b and s quark, respectively. We neglect the first term in 
eq. (^) since ^ 1 and p'^e^ = 0, which follows from ej^p^ = in the rest frame of the b 
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quark and the transversality condition e^fcf = 0, where is the transversal polarization vector 
of the ip meson f^. The second term can be written as 

— + 75)6 = — + 'J5)kif,b - isa^akiil + 75)&} ■ (H) 
ruh rrih 



Only the term in eq. ([Tl|) couples to the transversal component of the ijj and we obtain the 
corresponding amplitude as 

^(5 ^ s^T^ = -2CU{ml)^sa,^k'^Rbel, (12) 

^ mi) 

where R = ^-^^ denotes the chiral right projection. Note that the coupling structure is the 
same as due to a direct use of O7 = j^scTf^^mhRbF^'^ fl]] with the photon field strength tensor 
F'^'^ and rris = 0. For the tp"^ 7 conversion following the VMD mechanism we have 

< 0\J^,ei\i^^{h, e^)0 >= eQJ40)m^el , (13) 

where Qc = 2/3 and /^.(O) is the coupling at kf = 0, see eq. (|1^. Using the intermediate 
propagator of the ip meson at /c J = 0, we get 

A{b ^ si^^ ^ S7) = 2Cf40)^sa,^k'^Rbel . (14) 

The expression for the amplitude eq. ([1^) can be completed by summing over all cc resonant 
states V(15),^/'(25),^(3770),^/'(4040), ^^(4160) and ^(4415) 

A{b sijf ^ S7) = 2CY.fl{0)^sa,^ktRbel . (15) 

The various decay couplings = ftpi{jn%i) calculated using 

and the measured widths from [|1^] and given in Table ^ 

Here we need to extrapolate the coupling f^Xki = "^^'J /V'i(O)- We used the suppression 
factor [g 

'^ = 4(i5)(0)/4(i5)K) = 0.12 (17) 





Ih \GeV] 




U.4UO 


fi>{2S) 


n 989 


fi>{3770) 


0.099 


/V'(4040) 


0.175 


/i/>(4160) 


0.180 


/i/>(4415) 


0.145 



Table 2: Vector meson coupling constants used in the numerical calculations. 

obtained from data on the photoproduction of the and taking the factor k as universal for 
the other resonances. Q We now use eq. (|T5|) to find the matrix element of Bs 07 through 
the b — >■ sip'^ s'j transition at quark level. The matrix element is given as 



<^ip')\sa^aRk'^b\Bsip)> 



and we get the amplitude 



07) 



{4p.k,-p,h.ef)G{kl] 



2Cttt^^ 'J^eQS^t.up.k{p" 



rrib 



+ 9^i, 



m 



Bs 



-m(o), 



(19) 



where ei^, ef, are the polarization vectors and /ci, p' are the momenta of the photon and 
meson, respectively. We used Gikl = 0) = Fiikl = 0) [^]. Note, that the form factors 
introduced above are in general functions of two variables kf and p'^. We abbreviated here 
Fi{kl) = Fi{klp'^ = ml) and took Fi(0) = 0.24 ± 0.02 from g. 

Now we want to compare our result for A{Bs —>■ 07) eq. (19) with the same amplitude 
calculated by the method worked out in 0. This method is based on the new effective quark- 
gluon operator obtained by the interaction of the virtual charm quark loop with soft gluons, in 
contrast to a phenomenological description in terms of i{j resonances converting into a photon, 
as we used. In this approach, the operator Oi does not give any contribution to the matrix 
element of Bg —>■ 07 for an on-shell photon. The Fierz transformation of the operator O2 reads 



02 = l/NcOi + l/2 

octet ; 

^This is consistent with k = 0.11 |11| based on a dispersion relation calculatic 
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where 

Ooctet = 4:{c-r^^—YC){s-ff,^—Yb), (21) 

and A°/2 are the SU{3) colour generators. Then the only contribution comes from the colour 
octet part Ooctet- Using the operator Ooctet as a vertex of the virtual charm quark loop, which 
emits a real photon, and taking into account the c-quark-soft gluon interaction, a new effective 
operator is obtained. The matrix element of this operator between Bs and meson states gives 
the long distance amplitude of Bg 07 decay due to the O12 operators and it is written as 
(see for details; there the amplitude for the decay B K*^ is given) 

2 2 

A'{B, ^ 07) = 2C'e'te'^''{ze,,pXiP"' L + "^^^^^^g^^ L} , (22) 

where C = §^^2 '^g^ ■ The form factors L and L are calculated using QCD sum rules [0, 
rrib ml^ ml ^ 



,mn < ss > m? Airas < ss >^ ttI}, ml , m?,, 
[ Y2 27 + M^)]"^^^-M^)^ ' 

{7^{^ <G'> l„ds e-^ [-A + , ,,, (1 - -rfe)(l + -r^)] 



48 ' TT 'S M2+M'2' sM2 ' ' sM2 ' 

^ m2 <ss>ml 167ra^<_£s>^ 
- [ ]^ ^ + j^))- I2dj 

The Borel parameters M and M' are varied to find the stability region for L and L. We used 
the following input parameters in the sum rules: rn^ = 0.8 GeV^, < ss >= —0.011 GeV^, ^ < 
G2 >= 0.012 GeV\ = 0.23 GeV, fe, = 0.2 GeV, = 1.019 GeV and niB, = 5.369 GeV. 

The stabihty region is reached for 6 GeV^ < < 9 GeV^ and 2 GeV^ < M'^ < A GeV^ 
and we get 

L = (0.30 ± 0.05) GeV^ , 

L = (0.35 ±0.05) GeV^ . (24) 



Writing the amplitude for Bs 07 as 

A(')(B, ^ 07) = e^,e^^'{^e,,,X,p"^ ^"0 + g,, A+O) , (25) 



and using eq. (|T9D, ( ^21) and (^4]) , we can compare the coefficients obtained by two different 
methods and get 

\A- -A'-\ 



A- 
\A+-A'+\ 

This means, that the amphtudes agree within 10%. 



10% , 

5% . (26) 



In our approach, the structure for the transition b — > sip'^ 57 is proportional to ^-^kf 
(see eq. |12|), since the longitudinal part of the ip meson is disregarded to convert into a photon. 
Further, the form factors Fi{kf) and G{kf) in eq. ([T8|) are related for a real photon {kf = 0), 
Fi{0) = G{0). Therefore, in the amplitude A{Bs 07) appears only one form factor, which 
is Fi{0) in eq. (plQl). However, the form factors L and L in A'{Bs —>■ 07) given in eq. ( p2D are 
not related. They are calculated separately using QCD sum rules and this causes the difference 
between the ratios in eq. (|26|). In spite of the fact that the amplitudes A^ and A'^ are different 
from each other, they coincide within the given approximation and theoretical uncertainties 
lying in both methods. 

We can now present the amplitude for Bg 77 due to the chain reaction Bg — > (ptp — > 07 — >• 
77 using the intermediate propagator at zero momentum transfer and the — > 7 conversion 
vertex from the VMD model, 

< O|J^ez|0(p',e'^) >= eQsUiO)m^ef„ (27) 

where the polarization vector is treated as transversal. To apply the VMD mechanism to 
the amplitude eq. (|1^), we have to know the form factor at Fi{kl = 0,p'^ = 0). We took the 
extrapolated value -Fi(O) = -Fi(0, 0) = 0.16 ± 0.02 from p. Then the amplitude can be written 
with p' k2, — > €2 as 



A{Bs ^ 0V; ^ 07 ^ 77) = e'l{ki)e^2ik2)[g,.Al-^^^ + le^^o.^k'^k^A^-^J (28 



with the CP-even A'^ and CP-odd A parts 
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C - ^-^ o m i, — m 

At-Do, = 4x-i^i(0)E/4(0)eQe^^, 

^l-Do, = 4x^i^i(0)E/4(0)eQ- (29) 

where C is defined in eq. and the conversion factor x is defined as x = ~^Qs^^^- Here 



/^(O) = 0.18 GeV |TT| and Qs = —1/3. The extra factor 2 comes from the addition of the 



diagram with interchanged photons. 

The decay width for Bg 77 is obtained by adding the short distance amphtudes, the O7 
type LD effects and the LD effects due to the Bs (pip — > 07 — 77 transition, resulting in 

r(5, ^77) = — ^ — + +At-n p 

+ '^\A- + AIdo,+KdoS)- (30) 

3 Numerical estimates and discussion 

We presented a VMD model based calculation of the LD contribution to CP-even A'^ and 

CP-odd A^ decay amplitudes for Bg — > 77 decay due to the inclusive process h ^ sip. The 

conversions to photons from both the ipi resonances and the meson lead to two suppressions 

and make the amplitudes in eq. smaller compared to the ones from the LD effect of the 

operator O7 0. 

We estimated the ratio 
,+{-) 



and found 

2% < p < 4% (32) 

while varying ^ < /i < 2m6 and allowing 02 ^ to have a theoretical error of 25% as stated in 
1^. The analytical expression of the "effective" coefficient Cj^^ {p) of the operator O7 can be 
found in 
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As a conclusion, we investigated the LD-contribution to the 77 decay resulting from 

intermediate ipi production and compared our result with the one obtained by the interaction of 
the virtual charm loop with soft gluons 0. We see that both amplitudes are in good agreement 
within the errors of the calculation. The new LD contribution resulting from the four-quark 
operators Oi and O2 is smaller compared to the one coming from the Bg ^ </)'-/ ^ 77 chain 
decay 0] and effects our old estimate for the branching ratio B{Bs — > '~i'~i)sD+LDoj more 
than 1%. 
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